#include <iostream>
#include "interpolation_polynomial.h"
#include <vector>
using namespace std;

int main()
{
  vector<double> x = {0,0,3,3,5,5,8,8,13,13};
  vector<double> y = {0,0,225,225,383,383,623,623,993,993};
  vector<double> x_= {75,75,77,77,80,80,74,74,72,72};
  Hermit  p(x,y,x_);
  vector<vector<double>> f = p.getf();
  for(int i = 0; i < x.size();i++)
    {
      cout<< f[i][i]<< " ";
      cout<< endl;
    }
      cout<< endl;
  return 0;
}
